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Field of the Invention 

The present invention relates to image processing and, particularly multispectral 
imaging, and more particularly, to sharpening of images in the visual spectrum, especially 
those generated by airborne or satellite platforms. 

20 Copending Applications 

This application is copending with U.S. Patent Application Ser. No. , 



atty. docket number BOEI-1-1 181, and U.S. Patent Application Ser. No. , atty. 

docket number BOEI-1-1 182, all filed and all are hereby incorporated by 

reference. 

25 Background of the Invention 

Various systems exist for producing remotely sensed data for agriculture and other 
applications where it is desired to detect changes over time in an image area. One such 
sensing system is the Systeme Pour 1' observation de la Terre (SPOT) satellite developed by 
the French Centre National d'Etudes Spatiales (CNES). The SPOT satellite produces red, 
30 green, blue, and near-infrared color bands to a resolution of approximately 10 meters. For 
many uses, such as agriculture and other change detection applications, data produced by the 
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SPOT satellite meets various radiometric and temporal frequency requirements. However, 
this data is expensive. For example, a 60 km x 60 km image costs approximately $2,000. 

Another remote sensing system is a Land Remote Sensing Satellite Program 
(LANDSAT-7) satellite operated by NASA. The LANDSAT-7 satellite produces a single 
5 frame/array of a 1 80 km x 1 80 km image with red, green, blue, and near-infrared color bands 
at a resolution of about 25 meters and a panchromatic band that includes green, red, and near- 
infrared wavelengths at a resolution of 12.5 meters. A LANDSAT-7 image costs roughly 
$500 with a resulting cost and coverage benefit of 36. The data from the LANDSAT series 
satellite is attractive because of its affordable cost. However, the images produced by the 
10 LANDSAT series satellites do not have the multispectral spatial resolution necessary to 
produce high value agricultural or change detection information products. 

In order to overcome the low multispectral spatial resolution of the LANDSAT series 
satellite images, sharpening techniques have been introduced for increasing the resolution of 
images produced from the raw data generated by the LANDSAT series satellite. However, 
15 the images produced by these sharpening techniques lose some of the radiometric accuracy 
of the raw data, thereby reducing the value of the information produced. 

Therefore, there exists an unmet need in the art for increasing resolution of low cost 
data from remote sensing systems while preserving radiometric accuracy for effective use in 
agriculture and other change detection information products. 

20 Summary of the Invention 

The present invention provides a system, method, and computer program product for 
increasing resolution of low cost data from remote sensing systems while preserving 
radiometric accuracy for more highly effective use in agriculture and other change detection 
information products. Radiometrics is the detection and measurement of radiant 

25 electromagnetic energy. 

Embodiments of the present invention sharpen bands of sensor data in the visual 
spectrum. Blue, green, red, near-infrared, and panchromatic (pan) bands of data are received. 
Data of the pan band is corrected based on the blue, green, red, and near-infrared bands of 
data, and the data of one or more of the green, red, and near-infrared bands is sharpened 

30 based on the corrected data of the pan band. This technique can be applied to other 
multispectral bands covered by a pan band for the radiation region of interest, such as short 
wavelength infrared (SWIR), medium wavelength infrared (MWIR), or long wavelength 
infrared. 

In one aspect of the invention, the sharpened data of the green, red, and near-infrared 
35 bands is combined with the received pan band of data and then displayed. 
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In another aspect of the invention, the data suitably is corrected by converting the 
received blue, green, red, near-infrared, and pan bands of data to power-format, sharpening 
the converted power- formatted data of the green, red, and near-infrared bands based on the 
received pan band of data, and correcting the power- formatted data of the pan band based on 
5 sharpened power values for the data of the green, red, and near-infrared bands. 

In a further aspect of the invention, the data of the blue, green, red, and near-infrared 
bands may be resized in order to match the resolution of the data of the pan band prior to 
converting it to the power-format. 

In yet another aspect of the invention, the corrected pan band power suitably is 
10 subtracted from original pan band power to form a new band that covers red edge of 
vegetation. 

In still yet another aspect of the invention, the blue, green, red, near-infrared, and pan 
bands of data are generated by one of an aircraft or satellite sensing system as desired. 

Brief Description of the Drawings 

15 The preferred and alternative embodiments of the present invention are described in 

detail below with reference to the following drawings. 

FIGURE 1 is a block diagram of an exemplary system formed in accordance with an 
embodiment of the present invention; 

FIGURE 2 is an example of raw multi-spectrum data produced by a component of the 
20 system shown in FIGURE 1 ; 

FIGURES 3 and 4 show a flow diagram of an exemplary process performed by the 
system shown in FIGURE 1; and 

FIGURES 5A and 5B illustrate data progression with respect to the flow diagram of 
FIGURES 3 and 4. 

25 Detailed Description of the Invention 

By way of overview, the present invention provides a system, method, and computer 
program product for pan sharpening remotely sensed data, such as without limitation, 
LANDSAT series satellite data. The present invention performs power corrections and 
sharpening algorithms to the remotely sensed data in order to preserve spectral and spatial 

30 characteristics and features within an area of interest. 

FIGURE 1 illustrates an exemplary system 100 that performs the pan sharpening of 
remotely sensed data. The system 100 includes a remote sensing system 104, a processor 
106, a display device 108, and memory 110. The remote sensing system 104 senses data that 
includes a blue band, a red band, a green band, a near-infrared (NIR) band, and a 

35 panchromatic (pan) band. The pan band includes wavelengths of light that include green, red, 
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and NIR. The sensed data is sent to the processor 106, which sharpens the received data and 
makes it available for display on the display device 108 as described below. It will be 
appreciated that the process performed by the system 100 can be applied to other sets of 
bands that have pan bands which span the set of bands to be sharpened. 
5 FIGURE 2 illustrates exemplary data 120 produced by the remote sensing system 104 

(FIGURE 1). In one embodiment of the present invention, the remote sensing system 104 is a 
LANDSAT-7 satellite that produces a pan band 124 at a resolution of 15 meters. The 
LANDS AT-7 also senses data in a blue band 126, a green band 128, a red band 130, and a 
NIR band 132 at a resolution of 30 meters. Each of the color bands 126-132 are defined by 

10 width values w i5 w 2 , w 3 , and w 4 respectively, and are separated by gaps gi, g2, and g3, 
respectively. Other satellites that may also be used are Lconos or EO-1 with an ALI sensor. 

FIGURE 3 illustrates a non-limiting, exemplary process 150 performed by the 
system 100 (FIGURE 1). At a block 154, the pan band power of sensed pan band data is 
corrected. The sensed data is suitably a plurality of data units arranged in an array of data 

15 units. Each unit in the array has a power value that is being corrected in the block 154. 
Correction of the pan band power is described in more detail below in FIGURE 4. At a 
block 156, sensed data of the red, green, and NIR bands are sharpened based on the corrected 
pan band power. The sensed data of the green, red, and NIR bands are arrays of data units 
having a plurality of units of data. Each unit of data in the arrays is being sharpened at the 

20 block 156. At a block 158, the sharpened data of the red, green, and NIR bands are combined 
with the blue band and pan band to form an image. At a block 160, the image of the 
combined bands is displayed on the display device 108. At a block 164, the corrected pan 
band power is subtracted from the original pan band power to form a new band. The new 
band falls on the red edge of the vegetation spectrum. Changes that occur in this band can 

25 indicate various stresses on the vegetation including water stress. 

FIGURE 4 illustrates a non-limiting, exemplary process 200 for correcting the pan 
band power performed at the block 154 (FIGURE 3). At a block 204, the arrays of the blue, 
red, green, and NIR bands are resized to the resolution of the pan band. For example, the 
LANDSAT-7 satellite produces a data array of the pan band at a resolution of 15 meters, and 

30 data arrays of the blue, green, red, and NIR bands at 30 meters. The sensed data in each unit 
in the arrays are in digital number format. Each digital number of a unit for the blue, green, 
red, and NIR bands is divided into four subunits with each subunit having the same digital 
number value as the original unit. Other interpolation processes can be used to perform this 
rescaling process. As a result, each of the subunits is now at a resolution of 15 meters. At a 

35 block 206, the digital number of values of the resized data for the blue, green, red, and NIR 
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bands, and the data for the pan band are converted into radiance values. EQUATION (1) 
below shows an example formula for generating radiance values using calibration 
coefficients. 



y (LMAX T - LMAX T ) , _ , . x rj „ Ar 

5 T= - - — x (qcal - qcal mm) + LMIN T ( 1 ) 

(qcal max- qcal min) 

where 

Lt = spectral radiance in watts/(meter 2 * sr * micron) 

QCAL = quantized calibrated frame unit value in digital number (DN) format 
10 QCALMIN = minimized quantized calibrated frame unit value in DN format 

QCALMAX = maximum quantized calibrated frame unit value in DN format 
LMIN T = spectral radiance that is scaled to QCALMIN in watts/(meter 2 * sr * 
micron) 

LMAXt = spectral radiance that is scaled to QCALMIN in watts^meter 2 * sr * 
15 micron) 

At a block 208, each of the radiance values are converted into power values. 
EQUATION (2) below illustrates an exemplary formula for generating a power value from 
the radiance value. 

20 

P = Rad. • (wu - wi) (2) 
where 

w u = Upper width limit of associated band; and 
wi = Lower width limit of associated band. 
25 Rad = Radiance in band (this corresponds to L T ) 

The difference of the upper and lower width limits w u and wi for a band is multiplied 
by the radiance value to generate the power value. The upper and lower width limits w u and 
wi are determined by analysis of the spectral response of the respective band. 

30 At a block 210, the power values for each resized unit (subunit) of the red, green, and 

NIR band are sharpened. EQUATION (3) below illustrates an example of a power 
sharpening formula that determines a fraction of a power value P n in a unit in a band (green, 
red, or NIR) relative to the sum of the power of all the bands (green, red, and NIR) 
Pgreen + Pred + Pnir and multiplies the fraction by the original pan band power to generate a 

35 sharpened power value P s for the unit. 
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p 

P s = Original Pan Band Power x 2 (3) 

Pgreen + Pred + PnIR 

For example, a power value for one of the units for the green band is divided by the 
sum of the power values in the corresponding unit for the green, red, and NIR bands. This 
5 fraction is multiplied by the original pan band power value for the corresponding unit in 
order to produce a sharpened power value P s for the unit in the green band. 

At a block 212, width (w) and gap (g) information is derived from calibration data 
received from the remote sensing system 104. Referring back to FIGURE 2, width (w) and 
gap (g) information is shown relative to the green, blue, red, and NIR bands. At a block 214, 
10 the power of the pan band (i.e. power value in each unit of the pan band) is corrected using 
the power value for the corresponding unit in blue band, the sharpened power values for the 
corresponding units in the red, green, and NIR bands, and the derived width (w) and gap (g) 
information. EQUATION (4) below illustrates an exemplary formula for correcting the pan 
band power. 



15 



Pan Corrected Power = 



p pan 
(p2 + p3 + p4)x 2 



(p2 + p3 + p4)x 2 - p3 x f *£L + ^ 

^ w3v w3v ) 

_p4xf-^l-p2f-^+^]-plx glV 



w4v 



Vw2v w2v ) wlv 



(4) 

Where 

20 pi = multispectral power in band 1 (blue); 

p2 = multispectral power in band 2 (green); 

p3 = multispectral power in band 3 (red); 

p4 = multispectral power in band 4 (NIR); 

p pan = multispectral power in band 8 (panchromatic); 
25 wlv = difference between the wavelength lower limit for multispectral band 1 and 

the wavelength upper limit for multispectral band 1 ; 

w2v = difference between the wavelength lower limit for multispectral band 2 and 
the wavelength upper limit for multispectral band 2; 

w3v = difference between the wavelength lower limit for multispectral band 3 and 
30 the wavelength upper limit for multispectral band 3; 
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w4v = difference between the wavelength lower limit for multispectral band 4 and 
the wavelength upper limit for multispectral band 4; 

gl v = difference between the wavelength lower limit for multispectral band 2 and the 
wavelength upper limit for multispectral band 1 ; 
5 g2v = difference between the wavelength lower limit for multispectral band 3 and the 

wavelength upper limit for multispectral band 2; and 

g3v = difference between the wavelength lower limit for multispectral band 4 and the 
wavelength upper limit for multispectral band 3. 



10 The pan correction equation (4) subtracts out average power values of the gaps from 

the total pan power. 

If desired, the sharpened red, green, and NIR bands that are generated at the 
block 156 (FIGURE 3) can be reinserted into the process 200 in order to further improve the 
pan band corrected power. The processes 150 and 200 proceed until all units in a sensed 
15 array have been processed. It will be appreciated that reinserting the sharpened values into 
the process 200 is an interative process for refining the data and can be repeated as often as 
necessary. 

FIGURES 5A and 5B illustrate how the respective band data progresses through each 
of the steps of the processes 150 and 200. The process blocks are numbered the same as their 
20 corresponding blocks from FIGURES 3 and 4. Accordingly and for the sake of brevity and 
clarity, details of these process blocks need not be repeated for an understanding of the 
present invention. 

It will be appreciated that the present invention can operate in real-time, thereby 
producing sharpened images at video speed. 
25 While the preferred embodiment of the invention has been illustrated and described, 

as noted above, many changes can be made without departing from the spirit and scope of the 
invention. Accordingly, the scope of the invention is not limited by the disclosure of the 
preferred embodiment. Instead, the invention should be determined entirely by reference to 
the claims that follow. 
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